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Abstract. We study Gaussian valence bond states of continuous variable systems, obtained 
as the outputs of projection operations from an ancillary space of M infinitely entangled 
bonds connecting neighboring sites, applied at each of N sites of an harmonic chain. The 
entanglement distribution in Gaussian valence bond states can be controlled by varying the 
input amount of entanglement engineered in a (2M + l)-mode Gaussian state known as 
the building block, which is isomorphic to the projector applied at a given site. We show 
how this mechanism can be interpreted in terms of multiple entanglement swapping from 
the chain of ancillary bonds, through the building blocks. We provide optical schemes to 
produce bisymmetric three-mode Gaussian building blocks (which correspond to a single 
bond, M = 1), and study the entanglement structure in the output Gaussian valence bond 
states. The usefulness of such states for quantum communication protocols with continuous 
variables, like telecloning and teleportation networks, is finally discussed. 
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1. Introduction 

Quantum information aims at the treatment and transport of information in ways forbidden by 
classical physics. For this goal, continuous variables (CV) of atoms and light have emerged 
as a powerful tool (TJ. In this context, entanglement is an essential resource. Recently, the 
valence bond formalism, originally developed for spin systems [2], has been generalized to 
the CV scenario |f3j |4) for the special class of Gaussian states, which play a central role in 
theoretical and practical CV quantum information and communication (5). 

In this work we analyze feasible implementations of Gaussian valence bond states 
(GVBS) for quantum communication between many users in a CV setting, as enabled by 
their peculiar structure of distributed entanglement [4|. After recalling the necessary notation 
(Sec. I2J) and the construction of Gaussian valence bond states O (Sec. [3]), we discuss 
the characterization of entanglement and its distribution in such states as regulated by the 
entanglement properties of simpler states involved in the valence bond construction [4| 
(Sec. 3J. We then focus on the realization of GVBS by means of quantum optics, provide 
a scheme for their state engineering (Sec.|5]l, and discuss the applications of such resources in 
the context of CV telecloning |6, 7| on multimode harmonic rings (Sec. |6]>. 

2. Continuous variable systems and Gaussian states 

A CV system |fl~l[5) is described by a Hilbert space H = ^>^ =1 Hi resulting from the tensor 
product of infinite dimensional Fock spaces Hi's. Let and a\ be the annihilation and 
creation operators acting on Hi (ladder operators), and fa = (a^ + aj) and pi = (a^ — a\)/i 
be the related quadrature phase operators. Let R = (xi,pi, . . . , q^,p^) denote the vector of 
the operators fa and pi. The canonical commutation relations for the R4 can be expressed in 
terms of the symplectic form ft as 



The state of a C V system can be equivalently described by quasi-probability distributions 
defined on the 2 A-dimensional space associated to the quadratic form SI, known as quantum 
phase space. In the phase space picture, the tensor product H = Hi of the Hilbert spaces 
Hi's of the N modes results in the direct sum A = (J)^ Aj of the phase spaces Aj's. 

States with Gaussian quasi-probability distributions are referred to as Gaussian states. 
Such states are at the heart of information processing in CV systems [ 1 , 5 1 and are the subject 
of our analysis. By definition, a Gaussian state is completely characterized by the first and 
second statistical moments of the field operators, which will be denoted, respectively, by the 

vector of first moments R = (^(Ri), (R2), ■ ■ • , (R2N-1), (R2n)j an d me covariance matrix 
(CM) 7 of elements 



Coherent states, resulting from the application of displacement operators Dy = e lY nR 
(Y € R 2 ™) to the vacuum state, are Gaussian states with CM 7=1 and first statistical 
moments R = Y. First moments can be arbitrarily adjusted by local unitary operations 
(displacements), which cannot affect any property related to entropy or entanglement. They 



\Ri,Rj\ — 2i£lij 





7y = -(M-+4-^)-<^)(4-). 



(i) 
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can thus be assumed zero without any loss of generality. A TV-mode Gaussian state will be 
completely characterized by its real, symmetric, 2N x 2N CM 7. 

The canonical commutation relations and the positivity of the density matrix p of a 
Gaussian state imply the bona fide condition 

7 + in > , (2) 

as a necessary and sufficient constraint the matrix 7 has to fulfill to be a CM corresponding 
to a physical state |8, 9|. Note that the previous condition is necessary for the CM of any 
(generally non Gaussian) state, as it generalizes to many modes the Robertson-Schrodinger 
uncertainty relation 1 10 1. 

A major role in the theoretical and experimental manipulation of Gaussian states is 
played by unitary operations which preserve the Gaussian character of the states on which 
they act. Such operations are all those generated by terms of the first and second order in the 
field operators. As a consequence of the Stone-Von Neumann theorem, any such operation 
at the Hilbert space level corresponds, in phase space, to a symplectic transformation, i.e. to 
a linear transformation S which preserves the symplectic form n, so that n = S T nS, i.e. it 
preserves the commutators between the different operators. Symplectic transformations on a 
2iV-dimensional phase space form the (real) symplectic group, denoted by Sp( 2 N,xi)- Such 
transformations act linearly on first moments and "by congruence" on the CM {i.e. so that 
7 1—* SjS T ). One has Dot S = 1, V5 G Sp( 2 N,iR.)- A crucial symplectic operation is 
the one achieving the normal mode decomposition. Due to Williamson theorem [11 J, any 
A-mode Gaussian state can be symplectically diagonalized in phase space, so that its CM 
is brought in the form is, such that S~yS T = v, with v = diagj^i, v\,.,. fjv, vn}- The 
set {vi] of the positive-defined eigenvalues of \in~y\ constitutes the symplectic spectrum of 
7 and its elements, the so-called symplectic eigenvalues, must fulfill the conditions 14 > 1, 
following from the uncertainty principle Eq. and ensuring positivity of the density matrix 
p corresponding to 7. 

Ideal beam splitters, phase shifters and squeezers are described by symplectic 
transformations. In particular, a phase-free two-mode squeezing transformation, which 
corresponds to squeezing the first mode (say i) in one quadrature (say momentum, pi) and 
the second mode (say j) in the orthogonal quadrature (say position, qj) with the same degree 
of squeezing r, can be represented in phase space by the symplectic transformation 

Sij(r) — diagjexpr, exp— r, exp— r, expr} . (3) 

These trasformations occur for instance in parametric down conversions [12|. Another 
important example of symplectic operation is the ideal (phase-free) beam splitter, which acts 
on a pair of modes i and j as |fT3l 



B (9) / l_ * ^* C ° S ^ ^ S ^ n ^ 
y \ o,j 1— * en sin 9 — dj cos 9 



(4) 



and corresponds to a rotation in phase space of the form 

/ cos(0) shx(9) 

cos(0) sin(0) 

sin(0) -cos(0) 

V sin(6>) -cos(0) / 

The transmittivity r of the beam splitter is given by r = cos 2 (6) so that a 50:50 beam splitter 
(r = 1 /2) amounts to a phase-space rotation of ir/4. 
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o-ij (r) 



(6) 



The combined application of a two-mode squeezing and a 50:50 beam splitter realizes 
the entangling twin-beam transformation lfl4l 

T ij (r) = B ij (7r/4)-S ij (r), (5) 

which, if applied to two uncorrected vacuum modes i and j (whose initial CM is the identity 
matrix), results in the production of a pure two-mode squeezed Gaussian state with CM 

* iJ (r)=T ij (r)T?(r) given by 

/ cosh(2r) sinh(2r) \ 

cosh(2r) -sinh(2r) 

sinh(2r) cosh(2r) 

y -sinh(2r) cosh(2r) J 

The CV entanglement in the state (T\ j (r) increases unboundedly as a function of r, and 
in the limit r — > oo Eq. (|6]l approaches the (unnormalizable) Einstein-Podolski-Rosen (EPR) 
state 1151 , simultaneous eigenstate of relative position and total momentum of the two modes 
i and j. Concerning entanglement in general, the "positivity of partial transposition" (PPT) 
criterion states that a Gaussian CM 7 is separable (with respect to a 1 x N bipartition) if and 
only if the partially transposed CM 7 satisfies the uncertainty principle Eq. © (9] [16] • In 
phase space, partial transposition amounts to a mirror reflection of one quadrature associated 
to the single-mode partition. If {i>i} is the symplectic spectrum of the partially transposed 
CM 7, then a (N + l)-mode Gaussian state with CM 7 is separable if and only if P,* > 1 Vi. 
A proper measure of CV entanglement is the logarithmic negativity Ej^f ifTTl . which is readily 
computed in terms of the symplectic spectrum vi of 7 as 

Em = - Yj lo g^- (7) 

i: Vi<\ 

Such an entanglement monotone 1 18] quantifies the extent to which the PPT condition v. L > 1 
is violated. For 1 x JV Gaussian states, only the smallest symplectic eigenvalue £>_ of the 
partially transposed CM can be smaller than one [ 10], thus simplifying the expression of Ej^\ 
then the PPT criterion simply yields that 7 is entangled as soon as i>_ < 1, and infinite 
entanglement (accompanied by infinite energy in the state) is reached for — > + . 

For lxl Gaussian states 7^ j symmetric under mode permutations, the entanglement of 
formation Ep is computable as well via the formula |fl9l 

E F fr i j)=miac{0,f(i&)}, (8) 

with 

Being a monotonically decreasing function of the smallest symplectic eigenvalue v l _!_ J of the 
partial transpose j i .• of r ) i j, the entanglement of formation is completely equivalent to the 
logarithmic negativity in this case. For a two-mode state, Vi,j can be computed from the 
symplectic invariants of the state l20l . and experimentally estimated with measures of global 
and local purities 12~T1 (the purity \i = Tr p 2 of a Gaussian state p with CM 7 is equal to 
p= (Dct7)- 1 / 2 ). 



3. Gaussian valence bond states 



Let us review the basic definitions and notations for GVBS, as adopted in Ref. |4J. The 
so-called matrix product Gaussian states introduced in Ref. Q are JV-mode states obtained 
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Figure 1. Gaussian valence bond states. T tn is the state of N EPR bonds and -y is the three- 
mode building block. After the EPR measurements (depicted as curly brackets), the chain 
of modes -y x collapses into a Gaussian valence bond state with global state r° . See also 
Ref. (4). 



by taking a fixed number, M, of infinitely entangled ancillary bonds (EPR pairs) shared 
by adjacent sites, and applying an arbitrary 2M — + 1 Gaussian operation on each site 
i = 1, . . . , N. Such a construction, more properly definable as a "valence bond" picture 
for Gaussian states, can be better understood by resorting to the Jamiolkowski isomorphism 
between quantum operations and quantum states 11221 . In this framework, one starts with a 
chain of TV Gaussian states of 2M + 1 modes (the building blocks). The global Gaussian state 
of the chain is described by a CM T = ©^li 7^- As the interest in GVBS lies mainly in 
their connections with ground states of Hamiltonians invariant under translation 0, we can 
focus on pure (Det 7W = 1), translationally invariant (7W = 7 Vi) GVBS. Moreover, in this 
work we consider single-bonded GVBS, i.e. with M = 1. This is also physically motivated in 
view of experimental implementations of GVBS, as more than one EPR bond would result in a 
building block with five or more correlated modes, which appears technologically demanding. 

Under the considered prescriptions, the building block 7 is a pure Gaussian state of three 
modes. As we aim to construct a translationally invariant state, it is convenient to consider a 
7 whose first two modes, which will be combined with two identical halves of consecutive 
EPR bonds (see Fig. 01, have the same reduced CM. This yields a pure, three-mode Gaussian 
building block with the property of being bisymmetric [23 1, that is with a CM invariant under 
permutation of the first two modes. This choice of the building block is further justified 
by the fact that, among all pure three-mode Gaussian states, bisymmetric states maximize 
the genuine tripartite entanglement 11241 : no entanglement is thus wasted in the projection 
process. The 6x6 CM 7 of the building block can be written as follows in terms of 2 x 2 
submatrices, 



7 = 



Is 



sx 
£ sx 
Ix 



(9) 



The 4x4 CM of the first two modes (each of them having reduced CM 7 S ) will be denoted by 
7 SS , and will be regarded as the input port of the building block. On the other hand, the CM 
■y x of mode 3 will play the role of the output port. The intermodal correlations are encoded 
in the off-diagonal e matrices. Without loss of generality, we can assume 7 to be, up to local 
unitary operations, in the standard form [24] with 

7 S = diag{s, s}, ~i x = diagjz, x) , (10) 
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4^ 
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diag{w 



+ ■ 



-}; 



1 ± ^16s 4 - 8(x 2 + l)s 2 + (x 2 - If 



yj(x - 2s) 2 - 1 ± y/(x + 2s) 2 



The valence bond construction works as follows (see Fig. [3]). The global CM T — 
©fLi 7 acts as the projector from the state T m of the N ancillary EPR pairs, to the final 
TV-mode GVBS T out . This is realized by collapsing the state T ln , transposed in phase space, 
with the 'input port' T ss = ® . -f ss of T, so that the 'output port' T x = ® . ■y x turns into 
the desired T out . Here collapsing means that, at each site, the two two-mode states, each 
constituted by one mode (1 or 2) of 7 SS and one half of the EPR bond between site i and 
its neighbor (i — 1 or i + 1, respectively), undergo an "EPR measurement" i.e. are projected 
onto the infinitely entangled EPR state Il22l l3l. An EPR pair between modes i and j can be 
described, see Eq. ©, as a two-mode squeezed state <Tj,j(r) in the limit of infinite squeezing 
(r — ► oo). The input state is then T m = Hindoo (J)^ er^ where we have set periodic 

boundary conditions so that N + 1 = 1 in labeling the sites. The projection corresponds 
mathematically to taking a Schur complement (see Refs. [5] |22) for details), yielding an 
output pure GVBS of N modes on a ring with a CM 

-i r 



oy'Tsx, (ii) 

1, 1, —1} represents transposition in phase 



where T sx = 0^7^, and 9 = i diag{l, 
space (qi -> q i; pi -> -pi). 

Within the building block picture, the valence bond construction can be in toto 
understood as a multiple CV entanglement swapping ||251 . as shown in Fig. [3] the GVBS 
is created as the entanglement in the bonds is swapped to the chain of output modes via 
teleportation ll26l through the input port of the building blocks. It is thus clear that at a given 
initialization of the output port (i.e. at fixed x), changing the properties of the input port 
(i.e. varying s), which corresponds to implementing different Gaussian projections from the 
ancillary space to the physical one, will affect the structure and entanglement properties of 
the target GVBS. This link is explored in the following section. 



4. Entanglement distribution 

In Ref. [4 1 the quantum correlations of GVBS of the form Eq. (fTTT i have been studied, 
and related to the entanglement properties of the building block 7. Let us first recall the 
characterization of entanglement in the latter. As a consequence of the uncertainty principle 
Eq. (O, the CM Eq. (O of the building block describes a physical state if [24 1 

x + 1 

x > 1 , s > s min = - • (12) 

Let us keep the output parameter x fixed. Straightforward applications of the PPT separability 
conditions, and consequent calculations of the logarithmic negativity Eq. (0, reveal that the 
entanglement in the CM -f ss of the first two modes (input port) is monotonically increasing 
as a function of s, ranging from the case s = s m i n when 7 SS is separable to the limit s — > 00 
when the block 7 SS is infinitely entangled. Accordingly, the entanglement between each of 
the first two modes 7 S of 7 and the third one ~f x decreases with s. One can also show that the 
genuine tripartite entanglement in the building block increases with the difference s — s min 
ll24ll . The entanglement properties of the building block are summarized in Fig. [4] 
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Figure 2. How a Gaussian valence bond state is created via continuous-variable entanglement 
swapping. At each step, Alice attempts to teleport her mode (half of an EPR bond, depicted 
in yellow) to Bob, exploiting as an entangled resource two of the three modes of the building 
block (denoted at each step by 1 and 2). The curly bracket denotes homodyne detection, which 
together with classical communication and conditional displacement at Bob's side achieves 
teleportation. The state will be approximately recovered in mode 2, owned by Bob. Since 
mode 0, at each step, is entangled with the respective half of an EPR bond, the process swaps 
entanglement from the ancillary chain of the EPR bonds to the modes in the building block. 
The picture has to be followed column-wise. For ease of clarity, we depict the process as 
constituted by two sequences: in the first sequence [frames (1) to (4)] modes 1 and 2 are the 
two input modes of the building block (depicted in blue); in the second sequence [frames (5) 
to (8)] modes 1 and 2 are respectively an input and an output mode of the building block. 
As a result of the multiple entanglement swapping [frame (9)] the chain of the output modes 
(depicted in red), initially in a product state, is transformed into a translationally invariant 
Gaussian valence bond state, possessing in general multipartite entanglement among all the 
modes (depicted in magenta). 

The main question addressed in Ref. |4) is how the initial entanglement in the building 
block 7 redistributes in the Gaussian MPS r ou *. The answer is that the more entanglement 
one prepares in the input port j ss , the longer the range of pairwise quantum correlations in 
the output GVBS is, as pictorially shown in Fig. [4] 

In more detail, let us consider first a building block 7 with s = s m i„ = (x + l)/2. In 
this case, a separability analysis shows that, for an arbitrary number N of modes in the GVBS 
chain, the target state r ou * exhibits bipartite entanglement only between nearest neighbor 
modes, for any value of x > 1 (for x = 1 we trivially obtain a product state). In fact, each 
reduced two-mode block 7°"* is separable for |i — j\ > 1, 

With increasing s in the choice of the building block, one finds that in the target GVBS 
the correlations start to extend smoothly to distant modes. A series of thresholds s/. can be 
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Figure 3. Entanglement properties of the three-mode building block "y, Eq. (9), of the 
Gaussian valence bond construction, as functions of the standard form covariances x and 
d = s — s m j n . (a) Bipartite entanglement, as quantified by the logarithmic negativity, 
between the first two input-port modes 1 and 2; (b) Bipartite entanglement, as quantified by 
the logarithmic negativity, between each of the first two modes and the output-port mode 3; 
(c) Genuine tripartite entanglement, as quantified by the residual Gaussian contangle (2711241 . 
among all the three modes. 

found such that for s > Sk, two given modes i and j with \i — j\ < k are entangled. While 
trivially s\ (x) — s m ; n for any N (notice that nearest neighbors are entangled also for s = si), 
the entanglement boundaries for k > 1 are in general different functions of x, depending 
on the number of modes. We observe however a certain regularity in the process: Sk(x, N) 
always increases with the integer k. Very remarkably, this means that the maximum range 
of bipartite entanglement between two modes, or equivalently the maximum distribution of 
multipartite entanglement, in a GVBS on a translationally invariant ring, is monotonically 
related to the amount of entanglement in the reduced two-mode input port of the building 
block |4|. Moreover, no complete transfer of entanglement to more distant modes occurs: 
closer sites remain still entangled even when correlations between farther pairs arise. 

The most interesting feature is perhaps obtained when infinite entanglement is fed in 
the input port (s — > oo): in this limit, the output GVBS turns out to be a fully symmetric, 
permutation-invariant, TV-mode Gaussian state. This means that each individual mode is 
equally entangled with any other, no matter how distant they are [4]. These states, being thus 
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Figure 4. Pictorial representation of the entanglement between a probe (green) mode and its 
neighbor (magenta) modes on an harmonic ring with an underlying valence bond structure. 
As soon as the parameter s (encoding entanglement in the input port of the valence bond 
building block) is increased, pairwise entanglement between the probe mode and its farther 
and farther neighbors gradually appears in the corresponding output Gaussian valence bond 
states. By translational invariance, each mode exhibits the same entanglement structure with 
its respective neighbors. In the limit s —* oo, every single mode becomes equally entangled 
with every other single mode on the ring, independently of their relative distance: the Gaussian 
valence bond state is in this case fully symmetric. 



built by a symmetric distribution of infinite pairwise entanglement among multiple modes, 
achieve maximum genuine multiparty entanglement among all Gaussian states (at a given 
energy) while keeping the strongest possible bipartite one in any pair, a property known as 
monogamous but promiscuous entanglement sharing [27 1. 

Keeping Fig. [3] in mind, we can conclude that having the two input modes initially 
entangled in the building blocks, increases the efficiency of the entanglement-swapping 
mechanism, inducing correlations between distant modes on the GVBS chain, which enable to 
store and distribute joint information. In the asymptotic limit of an infinitely entangled input 
port of the building block, the entanglement range in the target GVBS states is engineered to 
be maximum, and communication between any two modes, independently of their distance, is 
enabled nonclassically. In the next sections, we investigate the possibility of producing GVBS 
with linear optics, and discuss with a specific example the usefulness of such resource states 
for multiparty C V quantum communication protocols such as telecloning [ 6 ] and teleportation 
networks lfl3ll . 

5. Optical implementation of Gaussian valence bond states 

The power of describing the production of GVBS in terms of physical states, the building 
blocks, rather than in terms of arbitrary non-unitary Gaussian maps, lies not only in the 
immediacy of the analytical treatment. From a practical point of view, the recipe of Fig.|3]can 
be directly implemented to produce GVBS experimentally in the domain of quantum optics. 
We first note that the EPR measurements are realized by the standard toolbox of a beamsplitter 
plus homodyne detection [22|, as demonstrated in several CV teleportation experiments |28|. 

The next ingredient to produce a TV-mode GVBS is constituted by N copies of the three- 
mode building block 7. We provide here an easy scheme (see also Refs. [6, 29]) to realize 
bisymmetric three-mode Gaussian states of the form Eq. (0. As shown in Fig. EJa), one 
can start from three vacuum modes and first apply a twin-beam operation to modes 1 and 3, 
characterized by a squeezing ri3, then apply another twin-beam operation to modes 1 and 2, 
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Figure 5. Optical production of bisymmetric three-mode Gaussian states, used as building 
blocks for the valence bond construction, (a) Three initial vacuum modes are entangled 
through two sequential twin-beam boxes, the first (parametrized by a squeezing degree ri3) 
acting on modes f and 3, and the second (parametrized by a squeezing degree r±2) acting on 
the transformed mode 1 and mode 2. The output is a pure three-mode Gaussian state whose 
covariance matrix is equivalent, up to local unitary operations, to the standard form given in 
Eq, {5). (b) Detail of the entangling twin-beam transformation. One input mode is squeezed 
in a quadrature, say momentum, of a degree r (this transformation is denoted by stretching 
arrows — <| |t— ); the other input mode is squeezed in the orthogonal quadrature, say position, of 
the same amount (this anti-squeezing transformation is denoted by the corresponding rotated 
symbol). Then the two squeezed modes are combined at a 50:50 beam-splitter. If the input 
modes are both in the vacuum state, the output is a pure two-mode squeezed Gaussian state, 
with entanglement proportional to the degree of squeezing r. 



parametrized by r±2 ■ The symplectic operation describing the twin-beam transformation (two- 
mode squeezing plus balanced beam splitter) is given by Eq. (0 and pictorially represented 
in Fig.|5jb). The output of this optical network is a pure, bisymmetric, three-mode Gaussian 
state with a CM 7s = T 12 {r 12 )T 13 (r 13 )T^(r 13 )T^{r 12 ) of the formEq. ©, with 



7. 



(e 4ri2 cosh (2r 13 ) + l) , -e~ 2ri2 (cosh (2r 13 ) + e 4ri2 ) 



diag { -e~ 2 ^ 



diag{cosh(2ri 3 ) , cosh(2ri 3 )} , 



(e 4ri2 cosh(2n 3 ) - l) , -e~ 2ri2 (cosh(2n 3 ) - e 4ri2 ) 



diag | V2e ri2 cosh (r i3 ) sinh (ri 3 ) , —V2e ri2 cosh (ri 3 ) sinh (ri 3 )j 



By means of local symplectic operations (unitary on the Hilbert space), like additional single- 
mode squeezings, the CM ~f B can be brought in the standard form of Eq. ( HOI , from which 
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For a given ri3 (i.e. at fixed x), the quantity r%2 is a monotonic function of the standard-form 
covariance s, so this squeezing parameter which enters in the production of the building block 
(see Fig.[5]l directly regulates the entanglement distribution in the target GVBS, as discussed 
in Sec. |4] 

The only unfeasible part of the scheme seems constituted by the ancillary EPR pairs. 
But are infinitely entangled bonds truly necessary? In Ref. [4| the possibility is considered of 
using a T ln given by the direct sum of two-mode squeezed states of Eq. ||6), but with finite 
r. Repeating the previous analysis to investigate the entanglement properties of the resulting 
GVBS with finitely entangled bonds, it is found that, at fixed (x, s), the entanglement in the 
various partitions is degraded as r decreases, as somehow expected. Crucially, this does not 
affect the connection between input entanglement and output correlation length. Numerical 
investigations show that, while the thresholds s& for the onset of entanglement between distant 
pairs are quantitatively modified - a bigger s is required at a given x to compensate the less 
entangled bonds - the overall structure stays untouched. This ensures that the possibility of 
engineering the entanglement structure in GVBS via the properties of the building block is 
robust against imperfect resources, definitely meaning that the presented scheme is feasible. 

Alternatively, one could from the beginning observe that the triples consisting of two 
projective measurements and one EPR pair can be replaced by a single projection onto the 
EPR state, applied at each site i between the input mode 2 of the building block and the 
consecutive input mode 1 of the building block of site i + 1 . The output of all the homody ne 
measurements will conditionally realize the target GVBS. 

6. Telecloning with Gaussian valence bond resources 

The protocol of CV quantum telecloning [6 1 among N parties is defined as a process in which 
one of them (Alice) owns an unknown coherent state, and wants to distribute her state to 
all the other N — 1 remote parties. The telecloning is achieved by a succession of standard 
two-party CV teleportations [26 1 between the sender Alice and each of the N — 1 remote 
receivers, exploiting each time the corresponding reduced two-mode state shared as resource 
by the selected pair of parties. The 1 — > 2 CV telecloning of unknown coherent states has 
been recently demonstrated experimentally JTJ . 

The no-cloning theorem [30] yields that the N — 1 remote clones can resemble the 
original input state only to a certain extent. The fidelity, which quantifies the success of a 
teleportation experiment, is defined as T = (ip m \p OMt \ip m ), where "in" and "out" denote the 
input and the output state. T reaches unity only for a perfect state transfer, p out = \i/j m )(ip' ln \ . 

Without using entanglement, by purely classical communication, an average fidelity of 
T c i = 1/2 is the best that can be achieved if the alphabet of input states includes all coherent 
states with even weight [311 . The sufficient fidelity criterion states that, if teleportation 
is performed with T > T c \, then the two parties exploited an entangled state [31]. The 
converse is generally false, i.e. some entangled resources may yield lower-than-classical 
fidelities. In Ref. 11321 it has been shown, however, that if the fidelity is optimized over all 
possible local unitary operations performed on the shared Gaussian resource (which preserve 
entanglement by definition), then it becomes equivalent, both qualitatively and quantitatively, 
to the entanglement in the resource. 



one has 
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Let us also recall that the fidelity of CV two-user teleportation [26 1 of arbitrary single- 
mode Gaussian states with CM -/ in (equal to the identity for coherent states) exploiting two- 
mode Gaussian resources with CM 7 b = ( £ab ) , can be computed ll33l as 

V £ ab lb J 

T = V5its ' S " 2 7m + ilai + lh + i£ab + ^ 1 (14) 

with £ = diag{ — 1 ,1}. We can now consider the general setting of 1 — > N — 1 telecloning, 
where N parties share a TV-mode GVBS as an entangled resource, and one of them plays the 
role of Alice (the sender) distributing imperfect copies of unknown coherent states to all the 
N — 1 receivers. For any N, the fidelity can be easily computed from the reduced two-mode 
CMs via Eq. ( TBI and will depend, for translationally invariant states, on the relative distance 
between the two considered modes. 

In this work we focus on a practical example of a GVBS on a translationally invariant 
harmonic ring, with N = 6 modes. As shown in the previous section, these states can be 
produced with the current optical technology. They are completely characterized, up to local 
unitary operations, by a 12 x 12 CM analytically obtained from Eq. ( fTTT i by considering the 
building block in standard form Eq. ©, whose elements are algebraic functions of s and x 
here omitted for brevity (as no particular insight is gained from their explicit expressions). 
First of all we can construct the reduced CMs 7°"+ fc of two modes with distance k, and 

evaluate for each k the respective symplectic eigenvalue D l _!_ l+k of the corresponding partial 
transpose. The entanglement condition s > Sk will correspond to the inequality y l ^ +k < 1. 
With this conditions one finds that S2 (x) is the only acceptable solution to the equation: 72s 8 — 
12(a; 2 + l)s 6 + (-34a; 4 + 28a; 2 -34)s 4 + (a; 6 -5a; 4 -5x 2 + l)s 2 + (a; 2 -l) 2 (a; 4 -6a; 2 + l) = 0, 
while for the next-next-nearest neighbors threshold one has simply S3 (x) = x. This enables us 
to classify the entanglement distribution and, more specifically, to observe the interaction scale 
in the GVBS T out : as discussed in Sec.|4]and explicitly shown in Ref. [4 |, by increasing initial 
entanglement in 7 SS one can gradually switch on pairwise quantum correlations between more 
and more distant sites. 

Accordingly, it is now interesting to test whether this entanglement is useful to achieve 
nonclassical telecloning towards distant receivers. In this specific instance, Alice will send 
two identical (approximate) clones to her nearest neighbors, two other identical clones (with 
in principle different fidelity than the previous case) to her next-nearest neighbors, and one 
final clone to the most distant site. The fidelities for the three transmissions can be computed 
from Eq. (TBI and are plotted in Fig. |6ja). For s = s m j n , obviously, only the two nearest 
neighbor clones can be teleported with nonclassical fidelity, as the reduced states of more 
distant pairs are separable. With increasing s also the state transfer to more distant sites is 
enabled with nonclassical efficiency, but not in the whole region of the space of parameters s 
and x in which the corresponding two-mode resources are entangled. 

As mentioned before, one can optimize the telecloning fidelity considering resources 
prepared in a different way but whose CM can be brought by local unitary operations (single- 
mode symplectic transformations) in the standard form of Eq. (fTTl i. For GVBS resources, 
this local-unitary freedom can be transferred to the preparation of the building block. A 
more general 7 locally equivalent to the standard form given in Eq. (fTOt . can be realized by 
complementing the presented state engineering scheme for the three-mode building block as in 
Eq. (foi l [see Fig. [3 a)], with additional single-mode rotations and squeezing transformations 
aimed at increasing the output fidelity in the target GVBS states, while keeping both the 
entanglement in the building block and consequently the entanglement in the final GVBS 
unchanged by definition. 
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Non— optimized fidelities 




Figure 6. 1 — * 5 quantum telecloning of unknown coherent states exploiting a six-mode 
translationally invariant Gaussian valence bond state as a shared resource. Alice owns mode 
i. Fidelities T for distributing clones to modes j such as k = \i — j\ are plotted for k = 1 
[(a),(d)]; k = 2 [(b), (e)]; and k = 3 [(c), (f)], as functions of the local invariants s and 
x of the building block. In the first row [(a)-(c)] the fidelities are achieved exploiting the 
non-optimized Gaussian valence bond resource in standard form. In the second row [(d)- 
(f)] fidelities optimized over local unitary operations on the resource are displayed, which 
are equivalent to the entanglement in the corresponding reduced two-mode states (see, as a 
comparison, Fig. 3 in Ref. (4)). Only nonclassical values of the fidelities (T > 0.5) are 
shown. 

The optimal telecloning fidelity, obtained in this way exploiting the results of Ref. [32|, 
is plotted in Fig.|6|b) for the three teleportations between modes i and j with k = \i — j\ = 
1, 2, 3. In this case, one immediately recovers a non-classical fidelity as soon as the 
separability condition s < s& is violated in the corresponding resources. Moreover, the 
optimal telecloning fidelity at a given k is itself a quantitative measure of the entanglement in 
the reduced two-mode resource, being equal to lF32l 

T° k pt = 1/(1 + !> M+fe ) , (15) 

where y" +k is the smallest symplectic eigenvalue of the partially transposed CM in 
the corresponding bipartition. The optimal fidelity is thus completely equivalent to the 
entanglement of formation Eq. ([8]) and to the logarithmic negativity Eq. (0. 

In the limit s — > oo, as discussed in Sec. 2J the GVBS become fully permutation- 
invariant for any N. Consequently, the (optimized and non-optimized) telecloning fidelity for 
distributing coherent states is equal for any pair of sender-receiverparties. These resources are 
thus useful for 1 — * N— 1 symmetric telecloning. However, due to the monogamy constraints 
on distribution of CV entanglement [27 1, this two-party fidelity will decrease with increasing 
N, vanishing in the limit TV — > oo where the resources become completely separable. In this 
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respect, it is worth pointing out that the fully symmetric GVBS resources are more useful 
for teleportation networks |fl3l [34], where N — 2 parties first perform local measurements 
(momentum detections) on their single-mode portion of the entangled resource to concentrate 
as much entanglement as possible onto the two-mode state of Alice and Bob, who can 
accomplish non-classical teleportation (after the outcomes of the N — 2 measurements are 
classically communicated to Bob). In this case, the optimal fidelity of TV-user teleportation 
network is an estimator of multipartite entanglement in the shared TV-mode resource [32|, 
which is indeed a GVBS obtained from an infinitely entangled building block. 

7. Conclusion 

The valence bond picture is a valuable framework to study the structure of correlations in 
quantum states of harmonic lattices. In fact, the motivation for such a formalism is quite 
different from the finite-dimensional case, where valence bond/matrix product states are 
useful to efficiently approximate ground states of TV-body systems - generally described by 
a number of parameters exponential in N - with polynomial resources [2]. In continuous 
variable systems, the key feature of GVBS lies in the understanding of their entanglement 
distribution as governed by the properties of simpler structures (4). This has also experimental 
implications giving a robust recipe to engineer correlations in many-body Gaussian states 
from feasible operations on the building blocks. We have provided a simple scheme to produce 
bisymmetric three-mode building blocks with linear optics, and discussed the subsequent 
implementation of the valence bond construction. We have also investigated the usefulness 
of such GVBS as resources for nonclassical communication, like telecloning of unknown 
coherent states to distant receivers on a harmonic ring. It would be interesting to employ the 
valence bond picture to describe quantum computation with continuous-variable cluster states 
ll35ll . and to devise efficient protocols for its optical implementation. 
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